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Diagnosing the Outflow from the SGR 1806—20 Giant Flare with 

Radio Observations 

J. Granot 1 , E. Ramirez-Ruiz 2 3 , G. B. Taylor 1,3 , D. Eichler 4 , Y. E. Lyubarsky 4 , 
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ABSTRACT 

On 2004 Dec. 27, the soft gamma repeater (SGR) 1806 — 20 emitted the 
brightest giant flare (GF) ever detected from an SGR. This burst of energy, 
which resulted in an (isotropic) energy release ~ 100 times greater than the only 
two other known SGR GFs, was followed by a very bright, fading radio afterglow. 
Extensive follow-up radio observations provided a wealth of information with un¬ 
precedented astrometric precision, revealing the temporal evolution of the source 
size, along with densely sampled light curves and spectra. Here we expand on 
our previous work on this source, by explaining these observations within one 
self-consistent dynamical model. In this scenario, the early radio emission is due 
to the outflow ejected during the GF energizing a thin shell surrounding a pre¬ 
existing cavity, where the observed steep temporal decay of the radio emission 
seen beginning on day 9 is attributed to the adiabatic cooling of the shocked 
shell. The shocked ejecta and external shell move outward together, driving a 
forward shock into the ambient medium, and are eventually decelerated by a 
reverse shock. As we show in Gelfand et al. (2005), the radio emission from 
the shocked external medium naturally peaks when significant deceleration oc¬ 
curs, and then decays relatively slowly. The dynamical modeling of the collision 
between the ejecta and the external shell together with the observed evolution 
of the source size (which is nicely reproduced in our model) suggest that most 
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of the energy in the outflow was in mildly relativistic material, with an initial 
expansion velocity v/c < 0.7di 5 , for a distance of 15di 5 kpc to SGR 1806 — 20. 
An initially highly relativistic outflow would not have produced a long coasting 
phase at a mildly relativistic expansion velocity, as was observed. 

Subject headings: pulsars: individual: (SGR 1806-20) - stars: neutron - stars: 
flare - stars: winds, outflows - hydrodynamics - ISM: bubbles 


1. Introduction 

Soft gamma repeaters (SGRs) are believed to be “magnetars” - a small class of slowly 
spinning neutron stars with extremely high surface magnetic fields, B > 10 15 G (Duncan & 
Thompson 1992; Kouveliotou et al. 1998). They have been discovered through their transient 
X-ray outbursts, during which they emit hundreds of short (~ 0.2 s), soft (kT ~ 25 keV) 
bursts; very rarely (only twice so far), SGRs emit giant flares (GF), extreme events with 
luminosities upwards of 10 44 erg s _1 . SGR 1806 — 20 lies in the Galactic plane, at a distance 
of about d = 15di5 kpc (Corbel & Eikenberry 2004; Cameron et al. 2005; McClure-Grifhths 
& Gaensler 2005). The third GF yet recorded occurred on 27 Dec. 2004, when SGR 1806 — 20 
emitted a burst so extreme that it was the brightest extra-solar transient event ever recorded 
(Palmer et al. 2005; Hurley et al. 2005). The flare was also unique in creating a very bright 
radio afterglow (Gaensler et al. 2005; Cameron & Kulkarni 2005), which was monitored for 
months, providing an amazing wealth of data at several radio frequencies, including the 
temporal evolution of the source size and shape, polarization and flux. 

The data from the radio source that appeared in the aftermath of the GF provide a rare 
opportunity for a detailed study of a mildly relativistic blast wave, which might help bridge 
the gap between the relativistic outflows in cosmological gamma-ray bursts and supernova 
remnants. Gelfand et al. (2005) present a rebrightening episode in the radio light curve, which 
is well fit by a semi-analytic spherical model where the radio emission resulted from a blast 
wave driven by > 10 24 ' 5 g of baryonic material driven off the neutron star, and the expanding 
radio nebula has now entered its Sedov-Taylor phase of evolution. An accompanying study 
of the evolution of the size of the nebula confirms that it is indeed decelerating (Taylor et 
al. 2005). Furthermore, the motion of the flux centroid implies a predominantly one-sided 
mildly collimated outflow, i.e. a wide one-sided jet (Taylor et al. 2005). 

In this paper, we expand on the framework laid out by Gaensler et al. (2005), Gelfand 
et al. (2005) and Taylor et al. (2005), and present a full dynamical model for the interaction 
of the outflow that was ejected during the 2004 Dec. 27 GF with its surroundings, and in 
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particular with an external shell. Our model, which is described in §2, explains the large 
and diverse data sets for this event and constrains the initial velocity of the ejecta from the 
flare. Both a relativistic (§2.1) and a Newtonian (§2.2) outflow are considered. We find that 
only a Newtonian outflow with an initial expansion velocity v ~ 0.7 d 15 c fits the observations 
well. In §3 we derive the synchrotron emission implied by our dynamical model and show 
that it also agrees nicely with the radio observations. Our conclusions are discussed in §4. 


2. The Underlying Dynamics 

The radio light curve initially exhibited a relatively moderate decay, ~ f -L5 , followed 
by an achromatic steepening at 4 ~ 9 days after the GF, to ~ t~ 2 '‘ (Gaenslcr et al. 2005). 
This was followed by a rebrightening starting at ~ 25 days and peaking at t p ~ 33 days 
(Gclfand et al. 2005); the decay rate slowed significantly thereafter (~ t L1 ; Gelfand et al. 
in preparation). 

The apparent expansion velocity, u ap = /3 ap c, of the nebula was initially fairly constant, 1 
~ 0.4di 5 c, and decreased at tdec ~ t p (Taylor et al. 2005). This value of v ap is for the 
geometrical mean of the semi-major axis and semi-minor axis, and assumes a double sided 
outflow. Under the latter assumption, u ap along the semi-major axis and semi-minor axis in 
~ 0.5di 5 c and ~ 0.25 di 5 c, respectively. However, as the motion of the flux centroid implies 
that the SGR is located close to the edge of the image along the semi-major axis (Taylor et 
al. 2005), this implies u ap ~ 1.0 d\ 5 c for the leading edge of the one-sided outflow. 

The true velocity (in units of c) of the emitting material near the edge of the image 
is f3 = (1 — W 2 ) 1 / 2 = /3 ap /(sind + fd^cosO) where 9 is the angle between its direction of 
motion and the line of sight. The minimal true velocity corresponding to a given observed 
value of (3 ap is /3 min = /3 ap ( 1 + /3 ap ) _1/2 or T min A™ = (3 ap and is obtained for an angle d min 
that satisfies cos 0 min = /3 min . The same angle, cos 9 = f3, is where /3 ap is maximal for a fixed 
(3 (and therefore /3 ap = Yj3 at the outer edge of the image of a spherical source expanding 
at a velocity (3 ). Therefore, we expect Yj3 ~ /3 ap , and in our case Yj3 ~ 1.0 d\ 5 , so that 
j3 ~ di 5 (l + d 2 5 )~ 1,/2 or (3 ~ 0.7 for di 5 « 1. 

Since the axis ratio of the radio image at the times relevant for our modeling is at most 
~ 2 : 1, and for the sake of simplicity, we adopt a spherical model for most of our analysis, 
and later introduce the corrections for a mildly collimated one sided outflow, as is implied 


1 We adopt the value that was derived by Taylor et al. (2005) for the average expansion velocity during 
the first 82 days, which is slightly higher than the value reported initially by Gaensler et al. (2005). 
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by the observations. 

As the source was already fading by the time of the first observation (ti ~ 7 days), 
the radio emission must have turned on at an earlier time and at a smaller radius. The 
radio spectrum and linear polarization observations suggest that we are seeing synchrotron 
emission. If the relativistic electrons that are emitting this radiation were accelerated at a 
much smaller radius (<C 10 16 cm) then most of their energy would have been lost via adi¬ 
abatic cooling by t\, thus dramatically increasing the GF energy requirements. Moreover, 
the achromatic light curve steepening at G (Gaensler et al. 2005), strongly suggests a hy¬ 
drodynamic transition at that time. A simple explanation for this overall behavior arises 
if the outflow from the GF initially coasted with negligible interaction with the ambient 
medium, until at t co i = 5 t co i i5 days it collided with a thin external shell, which caused it to 
decelerate by a reverse shock, while the external shell was accelerated by a forward shock. 2 
After this collision the two shells move together at a somewhat reduced speed. Thus the 
emission up to tb is dominated by the recently shocked electrons in these two shells. The 
radiation then arrives at the observer at a time t < 2 t co \ due to light travel effects and the 
finite time it takes for the shock to cross the shells. At t > tb the emission is dominated by 
the adiabatically cooling electrons in the two shells. As shown in §3 (see also Gaensler et al. 
2005) this naturally produces the unusually steep decay in the light curve. 

As the merged shell expands outwards it drives a shock into the ambient medium. An 
increasing amount of external mass is swept up, until the emission from this shocked fluid 
starts dominating the light curve at t > 25 days. This naturally produces a rebrightening in 
the light curve which peaks at t p ~ 33 days (Gelfand et al. 2005); as expected, a decrease in 
the expansion velocity was observed at about the same time, fd ec ~ t p (Taylor et al. 2005). At 
t > fdec the hydrodynamics gradually approach the self-similar Sedov-Taylor solution, which 
predicts a slower flux decay rate, in agreement with observations (Gelfand et al. 2005). An 
outline of our basic picture is shown in panel a of Fig. 1. Below we reproduce the main 
observed features using a simple analytic model for the interaction between the outflow and 
its surroundings, and present a numerical simulation that broadly agrees with the analytic 
model and nicely reproduces the observed evolution of the source size. 


2 Such a shell surrounding a pre-existing cavity is thought to be formed behind the bow shock due to 
the supersonic motion of SGR 1806 — 20 through the ISM and its quiescent wind (Gaensler et al. 2005). 
Alternatively, it could also arise from an earlier and initially faster mass ejection from the SGR 1806 — 20 GF, 
which was decelerated by the external medium to a velocity slightly below that of the coasting second shell, 
thus causing the two shells to collide (the “refreshed shock” scenario, e.g. Granot, Nakar & Piran 2003). 
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2.1. Relativistic Outflow 


A simple model for the collision between the cold ejecta shell of initial Lorentz factor 
r ej = (1 — /Jgj)^ 1 / 2 and mass M e j, and an external thin shell of mass M ext at rest at a radius 
_R ext , is a plastic collision where the two shells are shocked (the two shocked fluids separated 
by a contact discontinuity) and subsequently move together at Tj = (1 — Both T e j 

and T f are measured in the rest frame of the unperturbed external medium. Energy and 
momentum conservation in the rest frame of the merged shell require Ef/c 2 = Mf = Y r M e j + 
T /M ext and T r (3 r M ei = T f /3 f M ext , respectively, where E r = (1 - /3 2 )~ 1/2 = r ej r f (l - (3^(3 f) 
is the initial Lorentz factor of the ejecta in the rest frame of the merged shell. The resulting 
internal energy is E [nt /c 2 = (F r — 1 )M e j + (Tj — l)M ex t and the final velocity is 


Pf ( 1 1 ^ y 1 
Pej V r ej Mj 


( 1 ) 


This shows that an external mass of M ext ~ T e jM e j is required in order to significantly reduce 
the initial velocity. 

For an initially relativistic outflow T e j ~ E/M^c 2 S> 1 and /3 e j ~ 1 so that (3f ~ 
(1 + M ext c 2 /E)- 1 , which for M ext c 2 S> E (and correspondingly (3f <C 1) gives M ext vj ~ (3fE. 
Therefore, in this limit, the kinetic energy of the merged shell carries only a small fraction 
(~ (3f ) of the total energy, while most of the energy is in the internal energy of the shocked 
ejecta (E int ps E T e jM e jC 2 ). The relativistically hot shocked ejecta can then convert most 
of its internal energy back into kinetic energy through PdV work as the merged shell keeps 
expanding. This might initially (soon after the collision) accelerate the shell, and later cause 
it to decelerate more slowly with time (and radius), thus increasing the radius, Rd ec , where it 
decelerates significantly, compared to its value for a cold shell with the same (post collision) 
mass and velocity, Rd ec r^/ 2 1/3 Rexf 


Nevertheless, even if all the original energy is back in the form of kinetic energy at -Rdec, 
then still E /c 2 « M dec /3| ec « (3 f M ext where M dec = M(R dec ) « (47r/3)p ext i?^ ec is the total 
mass (in the shells and swept-up external medium) at -R dec , and /? dec = /?(-R dec ) ~ 0.4 d 15 > 
(3f. Finally, M ext ~ (47r/3)p ex t-Rext f° r m ost reasonable scenarios that produce an external 
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shell, such as a bow shock (Wilkin 1996). 3 This gives 


Rdec V Mdec ^ Pf 

Rex t / M ex t /3j ec 


C Q~ l 

~ r dec 


2 . 5 ^ 


( 2 ) 


We now proceed to compare the radio observations with the above calculations. The 
angular diameter of the source at the time of the first observation, t\ ~ 7 days, and at the 
epoch of deceleration, t dec ~ t p ~ 33 days, was about 80 mas and 300 mas, respectively. 4 The 
corresponding radii are Ri = 9.0 x 10 15 di5 cm and R dcc ~ 3.4 x 10 16 di5 cm. The requirement 
that R e .xt < R\ gives R dec /R ext > 3.75 and therefore (R dec /R ext ) 3 ~ 50(i?i/i? ex t) 3 > 50, 
which contradicts Eq. 2 for any reasonable value of dj 5 . Thus, an ultra-relativistic outflow 
(T e j 1) fails to reproduce the observations, since R dec would not be much larger than R ext ; 

_i jo 

specifically we expect R dec /R ext < lAd lb (see Eq. 2). 

2.2. Newtonian Outflow 

For a Newtonian outflow (/3 e j -C 1), Eq. 1 reduces to Pf/P e j ~ where Mf « 

M ej + M ext . Since M(R ext < r < R dec ) ~ M f and M ext « (47t/3 )p ext Rl xt , then M ext > M ej 
would imply R dec /R ext ~ 2 1/3 « 1.26, in contrast with observations. 

Therefore, we must have M e j S> M ext , which results in /5/ m /3 e j, Mf ~ M e j ~ 
M dec « (47r/3)pext-Rd ec , and M ei /M ext « (i? dec /i?ext) 3 > 50, or ifext « (t C ol/tl)i?I ~ 6.4 X 
10 15 f co i,5(ii5 cm and M e j/M ext « (7? de c/-Rext) 3 ~ 140t“ 3 5 . At f dec we directly measure the 
source size, i? dec , and expansion velocity, /3 e j ~ f3f ~ 7? d ec/ct de c ~ 0.4d 15 . Therefore, since 
the shocked external medium has comparable internal and kinetic energies, the energy in the 
outflow is given by E « (47T/3)p ex t ec ~ 3.8 x 10 46 nod'i 5 erg, and depends only on the un¬ 
known external density, n ext = p e xt/ m P = n 0 cm -3 . Here u e j = /3 e jc and E 46 = E/(10 46 erg). 


3 This is an important assumption. If somehow the mass of the shell would be larger by some factor 
/ = 3M ex t/47rp ex t-R 3 xt , then l?dec/-Rext would increase by a factor of / 1 / 3 so that / ~ 10 2 would be required 
in order to explain the observed evolution of the source size. Therefore, an external shell with / ~ 10 2 , 
or alternatively a sharp density drop by a factor of ~ 10 2 around R ex t that lasts for at least an order of 
magnitude in radius, would in principle be consistent with the observations. In practice, however, such 
external density profiles seem highly contrived and therefore not very likely. 

4 At both epochs the image is somewhat elongated and the quoted value is along the semi-major axis 
(Gaensler et al. 2005; Taylor et al. 2005). The ratio of the angular size at these two epochs, however, is 
rather robust and a comparable ratio is obtained along the semi-minor axis. This ratio is also applicable for 
a one sided outflow, as suggested in Taylor et al. (2005). 
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Thus n 0 ~ 0.26d^ 5 5 ^46 and M e j « E /r> 2 « 2.7 x 10 26 n 0 d? 5 g. These results for E and M ej 
are similar to those derived by Gelfand et al. (2005). 

A simple generalization for a wide one-sided jet is as follows. The volume of the shocked 
external fluid and therefore its mass for a given external density does not change. The kinetic 
energy per unit rest energy, T — 1, grows by a factor of ~ 4 — 5 at the head of the jet (where 
T/3 ~ 1, see beginning of §2) and decreases near the SGR. On average it increases by a factor 
of ~ 2 — 3, and the estimate for E/n ext increases by the same factor while M e Jn ext does not 
change. 

Finally, it is important to keep in mind that the outflow might consist of more than 
one component. The simplest example is a relativistic shell (with T > 1) followed by a 
Newtonian shell (with T/3 < 1) that was ejected slightly later during the GF. The relativistic 
shell is shocked and decelerated to a Newtonian velocity as it collides with the external 
shell, at t C ou, while the Newtonian shell catches up and collides with the slower merged 
relativistic + external shell at t co i ; 2 > Gold • As long as the velocity after the first collision is 
sufficiently smaller than that of the Newtonian shell, the subsequent dynamics would not be 
very different than for the Newtonian outffow case discussed above. An important difference, 
however, is that the emission would light up at ~ t co i,i/2r 2 <C Goi, 2 , he. much earlier than 
without the relativistic component. (A similar result is obtained also if there is a continuous 
external medium instead of a shell surrounding a cavity.) Rapid follow-up observations of 
future GFs could test this hypothesis directly, and teach us more about the properties of the 
outflow. In the present case, the later collision with the Newtonian shell might explain the 
change in the degree of linear polarization (from decreasing to increasing with time) and its 
position angle, at 7 10 days (Gaenslcr et al. 2005). 

We have tested the colliding shell scenario with the aid of numerical simulations, which 
model the dynamics much more accurately than the simple analytic model used above. Our 
basic picture is confirmed by these calculations, and the observed evolution of the source 
size is nicely reproduced (see Fig. 2). The simple analytic expression for the source size is 
FL(t > Goi) ~ Rdec niin[(t/t dec ), (t/t& e c ) 2//5 ]. The two asymptotic power laws correspond to the 
coasting phase and the Sedov-Taylor regime. A semi-analytic model which gives a smooth 
transition between these two asymptotic power laws is presented in Fig. 2 of Gelfand et al. 
(2005) and fitted to the source size as a function of time in Fig. 2 of Taylor et al. (2005). 
Around the time of the collision, t ~ 7 co i, there is a flat part in R(t) due to the shocks that 
are crossing both shells. This flat part is nicely reproduced by the numerical simulation, 
and its exact shape depends on the details of the collision (and thus does not have a simple 
and robust analytic description). In a future work (Ramirez-Ruiz et al., in preparation) we 
investigate the dynamics in more detail, including the implications of aspherical outflows, 



that are relevant given the elongated nature of the radio image (Gaensler et al. 2005) and 
the motion of its flux centroid (Taylor et al. 2005). 


3. Explaining the Observed Radio Emission 


Once the reverse shock crosses the ejecta and the forward shock crosses the external shell, 
the supply of newly accelerated electrons will be exhausted. As the merged shells expand, the 
existing relativistic electrons cool adiabatically and the magnetic held decreases, thus nicely 
reproducing the sharp decay that was observed in the radio flux between 9 and ~ 25 days, 
~ f -2 - 7 Gaensler et al. (2005). 

The emission from the forward shock is dominated by the newly shocked electrons 
which are accelerated to relativistic energies with a power law distribution, dN/d^ e oc , y~ p 
for 7 e > y m . At t < t dec the relative velocity of the shocked downstream fluid and the 
upstream fluid is roughly constant and equal to u re i = A-ei c ~ 0.3di 5 c since n re i/n e j ~ 3/4. 
The average Lorentz factor of the relativistic electrons is given by ( r j e ) = e e e' / ^en'm e c 2 where 
e' and n' are the proper internal energy density and number density of the shocked fluid, e e 
is the fraction of the post shock internal energy density in relativistic electrons, and £ e is the 
fraction of electrons that are accelerated to relativistic energies. The energy per proton is 
e' jn' = (T rel — 1 )m p c 2 ~ (/ 3 2 el /2)m p c 2 , where the second expression is valid in the limit of a 
Newtonian blast wave. For p > 2 we have y m = (7 e )(p — 2)/(p — 1) and therefore 


= £e / p- 2 \ rrip /3 2 el 
£ e \p - 1 ) m e 2 


2 g £ e ^e-l 



(3) 


where g = 3 (p — 2 )/{p — 1 ) (= 1 for p = 2.5), and e e _i = e e /0.1. Since 7 m is the lowest 
Lorentz factor of the relativistic electrons, by definition y m > 2. Gelfand et al. (2005) 
calculate the light curve under the assumption that e e /£ e = const and > 2 (see also Frail, 
Waxman & Kulkarni 2000), which is valid for e e > 0.1 or £ e -C 1 until there is significant 
deceleration. Once decreases to ~ 2 , the subsequent behavior of e e and depends on 
poorly understood shock acceleration of non-relativistic electrons. Here it is assumed that 
e e = const. Eq. 3 shows that for e ej _i < 1 (and it is difficult for e e to be much higher than 
0 . 1 ) we have y m ~ 2 which is constant all along, while £ e ~ (u re i/fej ) 2 (3 2 el decreases with 
time at t > t dec . This results in a more moderate temporal decay of the flux at t > f dec (see 
Eq. 4), in better agreement with observations. 

At f » tdec the shock dynamics approach the Sedov-Taylor self-similar solution, where 
R « (Et 2 /p ex t ) 1/5 . Therefore, n re i/n ej « (t dec /t)R/R dcc « min[l, (t/t d ec) _3/5 ] and v shfi /v ei « 
Vsh,o/vf = 4/3 where t dec = Rdec/^sh.o = (3E/27rp ex tV s 5 h,o) 1/3 - Here n sh ,o is the initial velocity 
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of the shock front for the blast wave propagating into the external medium. The post shock 
magnetic held is B = (Sire^e^t) 1 / 2 where e- mt = 2p ext v 2 el where e B = O.Ieb-i is the fraction 
of the post shock internal energy in the magnetic held. The number of synchrotron emitting 
electrons is N e = £ e M/m p where M = / 3) p ext R 3 and fb is the beaming factor (i.e. the 

fraction of the total solid angle occupied by the outflow). Finally, F u , max = N e P^ max /4nd 2 , 
where P u , max « P sy n/^syn, ^syn(7e) = (4:/3)a T c(B 2 / 8n)^, and is syn (%) = eB^ 2 e /2wm e c. The 
observed spectral slope in the radio suggest that we are in the spectral power law segment 
u m < v < u c , where u m = r / S y n (7m) and u c is the cooling break frequency. Thus we find 


F v = 4.2 f b g nQ (p+1)/20 P l |g 1+p)/10 e e> _ 1 


, 0 . 002 / 

at t > fdec, while F v [t < t dec ) « (t/t dec ) 3 Fv (Gee) • 


e B \ (p+ 1)/4 ( _ 2 ( V \ (l-p)/2 ( t 


d 


15 


(8.5 GHz) 


33 days 


-3(p+l)/10 


(4) 


The parameter values in Eq. 4 were chosen to match the observed flux at the peak of 
the rebrightening, Gee ps t p ~ 33 days. This demonstrates that an energy of r'U f b 10 46 erg, 
comparable to that in the GF (if it were emitted into a similar solid angle as the outflow), can 
be accommodated for reasonable values of the micro-physical parameters and the external 
density. Taking into account the relation n 0 ~ 0.26db5 5 -^46 derived in §2.2, we find that an 
equipartition limit e e ,e B < 0.3 (0.5) gives P 44 > 7.5 d 2 f (4.0 d 2 f), where P 44 = P/(10 44 erg), 
consistent with the conclusions of Gclfand et al. (2005). For a wide one sided jet of half¬ 
opening angle 9q ~ 0.5 rad, ft, = (1 — cos$o)/2 ~ 0.06, while E/n ext grows by a factor of 
~ 2 5 where E is the isotropic equivalent energy. Altogether, we obtain for the true energy 
P 44 > 5.7 d 2 f (3.4 d 2 f), as well as n 0 > 7.4 x lO” 3 ^ 2 ' 5 (4-4 x 10~ 3 d75 2 ' 5 ) and M e j j24 > 9.9 dff 
(5.9 dff), where M ej , 2 4 = M ej /(10 24 g). 


Finally, we estimate the expected flux at the end of the collision between the ejecta and 
the external shell, at ~ t co The external shell is accelerated to (3f ~ f3 e j while the ejecta are 
only slightly decelerated, so that the shock going into the external shell is stronger and likely 
to dominate the emission. The volume of the shell, 47r//i? 3 xt where rj = AR/R ext = O.lry-i, 
is reduced by a factor of 4 due to shock compression, and its internal energy is a fraction 
M ext /Mf ps Me X t/Mej ~ 0.007t 3 ol5 of the total energy E. Under similar assumptions as 
above, 


F„(t c oi) « 80 fbg p 7_i 


-1 „-(p+l)/4^(5+p)/4^,-l ^(p+l)/4^(5p—27)/4 


46 


e e-l e B,~l 


(1 


15 


V \ (l-U/2 


h5 GHz/ 


^coi.s m Jy , (5) 


in rough agreement with the extrapolation to t co i ~ 5 days of the observed flux, F v= 8.5 gh z = 
53 mJy, at the first epoch, t\ = 6.9 days (Cameron & Kulkarni 2005; Gaensler et al. 2005). 


For the parameter values used in Eq. 5 we obtain u m ~ 1 MHz, z/ sa ~ 50 MHz and 
v c ~ 10 1 ' Hz at t c oi ~ 5 days, where z/ sa is the self absorption frequency, so that the radio 


mJy 
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frequencies are well within the assumed power law segment of the spectrum. The low value 
we obtain for z/ sa is consistent with the lack of a change in the spectral slope down to 240 MHz 
(Cameron et al. 2005). Soon after the shock finishes crossing the shell, the electron power law 
energy distribution extends up to 7 max ~ 7 c (4oi) ~ 10 6 . Thereafter adiabatic cooling takes 
over and 7 max cx f -2 / 3 , while B oc f -1 so that zz syn (7max) ~ ^c(t co \)(t/t co i)^^ s . The emission 
from the shocked external medium starts to dominate at t ~ 25 days, i.e. t/t co \ ~ 5, and 
hence at that time r / syn (7max) > 10 15 Hz is well above the radio. 


4. Discussion 

We have described a dynamical model for the interaction with the surrounding medium 
of the outflow during the 2004 Dec. 27 giant flare (GF) from SGR 1806 — 20. This model 
nicely accounts for the observed radio light curves and spectrum as well as for the evolution 
of the source size with time. Using a simple analytic model, we have shown that the bulk 
of outflow from the GF could not have been highly relativistic, and was instead only mildly 
relativistic, with an average velocity of v < 0 . 7 di 5 C, similar to the observed roughly constant 
expansion velocity over the first month or so, taking into account the one-sided wide jet that 
is suggested by the radio data (Taylor et al. 2005). 5 6 

A major ingredient in our dynamical model is an external shell at a distance of 10 16 cm 
from the SGR. Such an external shell might naturally be formed by the bow shock due to the 
SGR’s quiescent wind and its supersonic motion relative to the external medium (Gaensler 
et al. 2005). A bow shock origin of the external shell has interesting implications. The bow 
shock stand-off radius is i?bs = 9-0 x lO 15 !/-^/^ cm, where u* = 250u 2 5o km s _1 is 

the velocity of SGR 1806 — 20 relative to the external medium, and L = 10 34 ' 5 L 34 .5 erg s -1 
is its spin-down luminosity . 6 In our scenario, 2i? ex t (the radius for a one sided outflow) is 
~ (1 — 2)i?bs, he. Rbs ~ (1 — 2)i?ext where R ext: ~ 6.4 x 10 15 t co i, 5 <ii 5 cm, which is the case for 


5 The local expansion velocity is highest along the jet axis. The observed axis ratio of ~ 2 : 1 for the 
radio image (Gaensler et al. 2005; Taylor et al. 2005) sets a lower limit on the true axis ratio of the emitting 
region. 

6 Before 1999 L « 8 x 10 33 ergs -1 while by 2001 and until before the Dec. 27th GF it leveled off at 
L ss 4.5 x 10 34 ergs -1 (Woods et al. in preparation). The dynamical time scale for the bow shock is 
tbs ~ Rbs/v * = 11 . 4 L 3 4 . 5 n 7 1 / 2 v 750 y r - I n our scenario, R hs ~ (1 - 2)i? ext so that t hs ~ I2t co b 5 di 5 v^ 0 yr. 
Since the spin down rate of SGR 1806 — 20 increased by a factor of ~ 5 several years before the GF, the 
steady state assumption for the bow shock is not valid for 1^250 <2 — 3. As a rough guide, one might still 
use the results for a steady wind (Wilkin 1996), with the average spin down luminosity over a period tbs- 
The exact shape of the bow shock could, however, be somewhat different than that of a steady wind. 
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our fiducial parameters. 

Lower limits on the energy, E > 10 44 ' 5 erg, and mass, M e j > 10 24 ' 5 g, of the outflow, 
and on the external density, n ext > 10 -2 cm -3 have been derived in §3 (see also Gelfand et 
ah 2005). The values of E, M e j and n ex t scale linearly with each other, E ~ 10 46 " 5 no erg 
and M e j ~ 10 26,5 no g. Note, however, that the minimal allowed density, no ~ 10~ 2 , requires 
u* ~ 2500 km s _1 and a similar kick velocity for SGR 1806 — 20 at its birth. While this is an 
extremely high kick velocity for typical radio pulsars, magnetars might have a significantly 
higher kick velocity compared to ordinary neutron stars, which might approach such high 
values (Duncan & Thompson 1992; Thompson & Duncan 1995). A lower kick velocity would 
suggest that the true values of E, M e j and n ext are larger than their lower limits by a factor 
of ~ l()0n 0 ~ 100AT 5 q. 

An alternative mechanism for producing the external shell is the ejection of a faster (and 
likely highly relativistic) component of the ejecta which carries a small fraction, /, of its total 
energy, just before (or in a slightly different direction relative to) the main part of the ejecta 
which carries most of its energy and is only mildly relativistic. The first component may be 
naturally identified with the matter that is coupled to the radiation in the initial spike of 
the GF, which is expected to reach a highly relativistic Lorentz factor (Nakar, Piran & Sari 
2005). In order to obtain the ratio Rd e J Rext ~ 4 — 5 that is inferred from observations, this 
would require / ~ (Rdec/Rext)~ 3 ~ 10“ 2 , i.e. that the ultra-relativistic component would 
carry ~ 1% of the total energy in the outflow. 

The fact that most of the energy in the outflow was in mildly relativistic ejecta, implies 
that the bulk of the ejecta was not coupled to the radiation of the initial spike. This could 
occur if the bulk of the outflow and of the radiation came out in different directions (either 
in small local patches, or within some global, possibly concentric, structure, as illustrated in 
panel b and c of Fig. 1, respectively). Thus, it is not even obvious whether they occupied a 
comparable solid angle, which is important when trying to compare their true energies. 

Our line of sight must have been in a relatively baryon-poor (and radiation rich) region, 
not only in order to see a bright initial spike, but also since otherwise the high optical depth 
of the electrons associated with the baryons in the outflow would have obscured the first 
~ 30 s of the GF tail emission. This can be seen as follows. If a shell of mildly relativistic 
proton-electron plasma with velocity vq and isotropic equivalent mass M e j, iso is ejected during 
the initial spike, at t ~ 0, then radiation emitted at time t after the initial spike would reach 
this shell at a radius R(t) = Vot/(l — Vq /c). The shell becomes optically thin to the radiation 
from the tail of the GF the a time f* when Tt(G) = Mej i i SO cr T /[m p 47ri? 2 (t*)] = 1, i.e. at 
f* = (1 — Vq /c) Vq 1 ( M e jd SO (jT /47rm p ) 1 A. This gives f* > 25 s from the lower limits we obtain 
for M ej , iso and for /3 min « 0.7 for a one-sided outflow (a similar result is obtained for the 
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spherical case). The presence of e ± pairs in the outflow or a large mass or velocity of the 
ejecta would only increase f*. 

In order for the highly relativistic ejecta to form the external shell into which the 
mildly relativistic ejecta collides, the two components must have a significant overlap in 
solid angle by the time of the collision, i.e. after the ultra-relativistic component decelerates 
to Newtonian velocities. This could occur if the two components form small patches, rather 
than a large scale coherent structure (such as a baryon-rich core surrounded by a baryon- 
poor and radiation-rich outer region). Alternatively, if the external shell was created by the 
bow shock, then an energy ratio of / < 1CU 2 is required in order for the highly relativistic 
component not to alter the radio emission considerably. This would be consistent with the 
very low energy in the ultra-relativistic component that is expected for a pure pair plasma 
(or for a very low baryon loading; Nakar, Piran & Sari 2005). 

For a wide one-sided jet we find that the lower limit on the isotropic equivalent ki¬ 
netic energy in the outflow, > 5 x 10 45 d 2 5 5 erg, is only a factor of rsj (3 — 4)d 15 0 ' 5 smaller 
than the isotropic equivalent energy radiated in the GF itself. This suggests that the two 
isotropic equivalent energies are comparable. If the solids angles occupied by the baryon-rich 
(radiation-poor) and by the radiation-rich (baryon-poor) regions are comparable (which is 
not at all obvious), then this would suggest that the true energies in the GF itself and in 
the kinetic energy of the outflow are of the same order. 

A much dimmer radio afterglow was detected following the 1998 Aug. 27 GF from 
SGR 1900 + 14 (Frail, Kulkarni & Bloom 1999), which despite the much sparser data, shows 
similarities to the radio afterglow discussed here. This suggests that our model might be 
applicable more generally. The spin down luminosity L of the two SGRs is comparable, and 
so is the time at which the light curve started to decay steeply (~ 9 — 10 days), suggesting 
that i?bs hr each case is not very different. This would imply a similar n ex t v 2 . Under these 
assumptions, the large difference in the radio luminosity (by a factor of ~ 500) would be 
mainly a result of the much larger energy content carried by the outflow of SGR 1806 — 20 
immediately after the GF. 
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Fig. 1.— An illustration of the basic underlying geometry in our model, a) A pre-existing 
shell surrounding a cavity (i.e. an evacuated region) is possibly formed due to the interaction 
of the SGR quiescent wind with the external medium, and the SGR’s super-sonic motion 
relative to the external medium. The outflow from the SGR 1806 — 20 Dec. 27 giant flare 
was ejected mainly in one preferred direction, probably not aligned with the head of the 
bow shock (which is in the direction of the SGR’s systemic motion). The ejecta collide 
with the external shell at a radius R ext and then the merged shell of shocked ejecta and 
shocked swept up external shell keeps moving outward at a constant (mildly relativistic) 
velocity. As it coasts outward, it gradually sweeps the external medium until at a radius 
Rdec ~ (4 — 5) Rext it sweeps up a sufficient mass of external medium in order to decelerate 
significantly. At R > Rd cc the structure of the flow gradually approaches the spherical self¬ 
similar Sedov-Taylor solution. b,c) Most of the mass in the outflow was in baryons that were 
decoupled from most of the radiation, and our line of sight was baryon poor. This naturally 
occurs if there are separate baryon rich (radiation poor) and baryon poor (radiation rich) 
regions. Such regions might consist of small baryon rich clumps surrounded by baryon poor 
regions (panel b), or might alternatively be part of a global large scale, possibly concentric, 
configuration (panel b). 
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Fig. 2.— Temporal evolution of the observed size of the source from radio measurements 
( asterices ) of Taylor et al. (2005) assuming d 15 = 1, together with the source size from 
our numerical simulation ( triangles ) and for our simple analytic model (dotted line). Our 
numerical simulation featured the collision between the outflow ejected during the SGR 
giant flare and a pre-existing shell surrounding a cavity. The calculations were done in two- 
dimensional cylindrical coordinates for ten levels of refinement using the piecewise parabolic 
method (PPM) adaptive mesh refinement code FLASH. The spherical initial configuration 
is as follows. In the inner region (outflow from the SGR, inner 5 x 10 14 cm) both a thermal 
energy of E = 10 46 erg and ejecta mass, M e j, are distributed uniformly; M e j is selected so 
that v = (2E/M e j) 1 / 2 ~ 0.4c (i.e. M e j ~ 1.4 x 10 26 g). The injected gas and surrounding 
ISM (with p ext = 2 x 10~ 24 g cm -3 ) are characterized by a 5/3 adiabatic index. More details 
will be presented in Ramirez-Ruiz et al. (2005, in preparation). 











